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ABSTRACT 

We report here on 5 GHz global very-long-baseline interferometry (VLBI) observations 
of SN 1986J, 16 yr after its explosion. We obtained a high-resolution image of the 
supernova, which shows a distorted shell of radio emission, indicative of a deformation 
of the shock front. The angular size of the shell is ~ 4.7 mas, corresponding to a linear 
size of ~ 6.8 x 10 17 cm for a distance of 9.6 Mpc to NGC 891. The average speed 
of the shell has decreased from ~7400 km s _1 in 1988.74 down to about 6300 km s -1 
in 1999.14, indicative of a mild deceleration in the expansion of SN1986J. Assuming 
a standard density profile for the progenitor wind (p cs ex r -s ,s = 2), the swept-up 
mass by the shock front is ^2.2 Mq. This large swept-up mass, coupled with the mild 
deceleration suffered by the supernova, suggests that the mass of the hydrogen-rich 
envelope ejected at explosion was ^ 12 M Q . Thus, the supernova progenitor must 
have kept intact most of its hydrogen-rich envelope by the time of explosion, which 
favours a single, massive star progenitor scenario. We find a flux density for SN 1986J 
of ^7.2 mjy at the observing frequency of 5 GHz, which results in a radio luminosity 
of - 1.4 x 10 37 ergs" 1 for the frequency range 10 7 -10 10 Hz (a = -0.69; S v cx v a ). We 
detect four bright knots that delineate the shell structure, and an absolute minimum 
of emission, which we tentatively identify with the centre of the supernova explosion. 
If this is the case, SN 1986J has then suffered an asymmetric expansion. We suggest 
that this asymmetry is due to the collision of the supernova ejecta with an anisotropic, 
clumpy (or filamentary) medium. 

Key words: Techniques: interferometric - supernovae: individual: SN 1986J - ISM: 
supernova remnants - Radio continuum: stars - Galaxies: individual: NGC 891 



1 INTRODUCTION 

Type II supernovae (SNe) are associated with massive stars 
that have expelled slow, dense winds during their super- 
giant phase. The stellar explosion drives a shock into this 
wind, at speeds as high as 20000 kms -1 and temperatures of 
~ 10 9 K. In addition, a reverse shock propagates back into 
the stellar envelope at speeds of 500-1000 kms -1 relative 
to the expanding ejecta. This is the so-called standard in- 
teraction model (hereafter SIM; Chevalier 1982, Nadyozhin 
1985), and radio, optical, and X-ray emission from Type II 
supernovae have been usually interpreted within this model. 
The outgoing shock forms a high-energy-density shell that 
is responsible for the production of synchrotron radio emis- 
sion, while the reverse shock accounts for the optical and 
soft X-ray emission. 



SN 1986J in NGC 891 is one of the most radio luminous 
SNe ever discovered. Indeed, at a distance of ~ 9.6 Mpc 
(Tully 1998), it had a peak luminosity at 5 GHz about 8 
and 13 times greater than SN 1979C and SN 1993J, respec- 
tively. The precise date of its explosion is not known, but on 
the basis of the available radio and optical data SN 1986J 
was estimated to have exploded around the end of 1982, or 
the beginning of 1983 (Rupen et al. 1987, Chevalier 1987, 
Weiler, Panagia & Sramek 1990). Based u pon it s large ra- 
dio luminosity, Weiler, Panagia, & Sramek (199C) suggested 
that the progenitor star was probably a red giant with a 
main-sequence mass of 20 - 30 Mq that had lost material 
rapidly ( M £ 2 x 10" 4 M yr ') in a dense stellar wind. 



VLBI observations by Bartel et al. (1991) showed that the 
radio structure of SN1986J had the form of a shell, in agree- 
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ment with expectations from the SIM model, with the min- 
imum of emission located approximately at its center. 

SN 1986J is peculiar in several respects, (i) Bartel et al. 



(1991) showed the presence of protrusions in the brightness 
distribution of the supernova. This finding was interpreted 
as evidence of deviation from spherical symmetry. Another 
supernova for which considerable asymmetry in its radio 
structure has been found is SN1987A (Gaensler et al. 1997, 
and references therein). Those cases contrast with the re- 
markable spherically symmetric expansion observed for the 
radio shell of SN 1993J (Marcaide et al. 1995, 1997; Bartel 
et al. 2000); (ii) the radio light curves for SN 1986 J were 
not well fitted within the SIM for radio supernovae, which 
led Weiler et al. (1990) to invoke the existence of a mixed 
medium (thermal absorbers and non-thermal emitters), or 
significant filamentation in the circumstellar medium; (iii) 
optical spectra obtained in 1986 (Rupen et al. 1987) and in 
1989 (Leibundgut et al. 1991) showed optical narrow lines 
with FWHM « 500 - 700 kms -1 . In the SIM, the optical and 
X-ray emission arise in the reverse shock, which moves with 
a speed close to that of the forward shock, i.e., several thou- 
sands kms" 1 . Since the optical velocities are far lower than 
those determined from VLBI (radio) observations at simi- 
lar epochs, these narrow lines question the validity of the 
SIM to successfully explain the ob serve d optical and X-ray 
emission from SN1986J. Chugai (1993) proposed a differ- 
ent model which overcomes the problems of the SIM and 
seems plausible for SN 1986J. In his model, the supernova 
envelope is not colliding with a smooth distribution of wind 
material, but with a clumpy one, and the bulk of the ob- 
served X-rays originates in the shocked dense wind clumps. 
Therefore, the narrow-line (Av ~ 500 kms -1 ) material is 
not directly related to the shock wave (which moves at much 
higher speeds), but to the speed of the s hock-e xcited dense 
clouds. More recently, Chugai & Belous (1999) have shown 



that the evolution of the radio light curves of SN1986J can 
be well interpreted by free- free absorption in a clumpy wind. 
The possible presence of clumps of radio absorbing plasma 
is also consistent with the proposal by Weiler et al. (1990) 
for a filamentary structure in the circumstellar medium. 

Unfortunately, VLBI observations of SN 1986J have 
been so scarce since its explosion that only one high- 
resolution radio image for SN 1986J exists (Bartel et al. 
1991). In this paper, we use very sensitive 5 GHz global 
VLBI observations taken on 21 February 1999 to obtain the 
second high-resolution radio image of SN 1986J, 16 yr after 
its explosion. We present a brief report of the VLBI observa- 
tions in section ^| present and discuss our results in section 
and summarize our conclusions in section ^. 



2 VLBI OBSERVATIONS AND IMAGE 
PROCESSING 

We used archival global VLBI data for SN1986J. The su- 
pernova was observed at a frequency of 5 GHz from 17:00 
UT on 21 February 1999 to 04:50 UT on 22 February 1999, 
using a very sensitive VLBI array that included the follow- 
ing antennas (diameter, location): The whole VLB A (25m, 
10 identical antennas across the US), phased- VLA (130m- 
equivalent, NM, US), Effelsberg (100m, Germany), Medicina 
and Noto (32m each, Italy), and Onsala (20m, Sweden). The 



telescopes received both left- and right-hand circular polar- 
izations (LCP and RCP) which, after correlation, were com- 
bined to obtain the total intensity mage presented in this pa- 
per. Effelsberg had technical problems (damaged gear) and 
thus could not take part in the observations. Onsala only 
recorded in LCP mode, and therefore its data were not used 
to obtain our total intensity VLBI image. 

The observations were made with a bandwidth of 64 
MHz. The data were correlated at the VLBA Correlator of 
the National Radio Astronomy Observatory (NRAO) in So- 
corro (NM, US). The correlator used a pre-averaging time 
of 3 s. Since SN 1986J was expected to be very faint, at 
a level of a few mjy, the observations were carried out 
in phase-reference mode. SN 1986J and the nearby ICRF 
source, 3C 66A, were alternately observed during the 12-hr 
experiment. The observations consisted of ~ 120 s scans on 
SN1986J and of ~ 70 s scans on 3C66A, plus a few addi- 
tional seconds of antenna slew time to make a duty cycle of 
190 s. 3C 66A was also observed as amplitude calibrator for 
SN 1986J. The source 3C 84 was observed as a fringe finder. 

The correlator data were analyzed using the Astronom- 
ical Image Processing System (AIPS). The visibility ampli- 
tudes were calibrated using the system temperature and gain 
information provided for each telescope. The instrumental 
phase and delay offsets among the 8-MHz baseband con- 
verters in each antenna were corrected using a phase cal- 
ibration determined from observations of 3C 84. The data 
for the calibrator 3C 66A were then fringe-fitted in a stan- 
dard manner. We exported the 3C 66A data into the Caltech 
imaging program DIFMAP (Shepherd et al. 1995) for map- 
ping purposes. 3C 66A showed a flux of ~ 0.95 Jy at 5 GHz, 
and displayed a one-sided core-jet structure at mas scales, 
with the jet extending southwards up to ~ 28 mas. The fi- 
nal source model obtained for 3C 66A was then included as 
an input model in a new fringe-fitting search for 3C 66A. 
In this way, the solutions obtained were structure-free. The 
phases, delays, and delay-rates determined for 3C 66 A were 
then interpolated and applied to the source SN 1986J. The 
SN 1986J data were then transferred into the DIFMAP pro- 
gram. Standard self-calibration techniques and purely uni- 
form weighting were used to achieve maximum resolution in 
the hybrid maps shown in Fig. jjj. 

We note two powerful aspects of the use of the phase- 
reference technique: (i) it effectively increases the integra- 
tion time on a source from minutes to hours, thus allow- 
ing the detection and imaging of very faint objects (e.g. 
Beasley and Conway 1995), and (ii) it retains the positional 
information of SN 1986J with respect to the phase-reference 
source, 3C66A, almost 1 degree apart. Indeed, a Fourier 
inversion of the SN 1986J data within DIFMAP, before ap- 
plying any self-calibration, showed that the a priori posi- 
tion of SN1986J used in the correlator (a = 02 1 ? 22^3 1=320, 
S = +42?19;57'.'282; J2000.0) was off from the actual one by 
~17 mas in right ascension (westwards), and by ~21 mas 
in declination (southwards). Therefore, phase-referencing 
proved to be also useful in providing further refinements 
in the precise coordinates of SN 1986 J. 
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Figure 1. (a) Hybrid map of the galaxy NGC 891 and its supernova SN 1986J made with the Very Large Array (VLA) at the same 
frequency (5 GHz) and epoch of the global VLBI observations on February 21, 1999. The contours are (3,5,10,20,40,80,160,300) X 25 
fi3y beam -1 , the root-mean-square (rms) noise off-source. The peak of brightness of the map corresponds to the supernova and is ~ 8.2 
mjy beam -1 . The dimensions (full width at half maximum) of the restoring beam are 7.4 X 4.0 arcsec 2 , with the beam's major axis 
oriented along a position angle of 88? (b) Global Very-Long-Baseline Interferometry (VLBI) hybrid map of SN 1986J on February 21, 
1999. The contours are (3,5,7,9,11,13,15,17,19) X 56 fijy beam -1 , the rms noise off-source. The peak of brightness of the map is 1.13 
mjy beam -1 and the restoring beam (bottom left in the image) is 1.3 X 0.9 mas 2 at a position angle of -13? 4. In both panels, north is 
up and east is left. 



3 RESULTS AND DISCUSSION 

The only other available radio image of SN 1986J is that ob- 
tained by Bartel et al. (1991; hereafter B91) of 29 September 
1988, at a frequency of 8.4 GHz. Our global VLBI image of 
21 February 1999 was obtained at 5 GHz. The most remark- 
able thing in our image (panel (b) in Fig. [j]) is the fact that 
the brightness distribution of SN 1986J still shows the form 
of a shell, even though more than 10 yr since the last obser- 
vations have passed. However, the shell, definitely distorted 
and with strong asymmetries, shows changes in the bright- 
ness distribution between the two epochs. 



3.1 Asymmetric shell morphology and protrusions 
in SN 1986J 

We distinguish at least four bright features in the rim of the 
shell labeled A,B, C, and D in Fig. |l|. The absolute mini- 
mum (~0.26 mjy beam -1 ) of the shell brightness distribu- 
tion is indicated by the X in panel (b). The flux densities 
and distances of the bright features A,B,C, and D from X 
are shown in Table [j]. At a flux density of 1.13 mjy beam -1 , 
the corresponding brightness temperature of feature A is 
- 1.8 X 10 7 K. Features B and C are local maxima with 
brightness temperatures of ~1.0 and ~1.2 x 10 7 K, respec- 
tively. Component D is an extended feature in the AC di- 



rection. Therefore, its flux density and radial distance from 
X are just indicative. 

The image by B91 showed several 'knots' along the rim 
of the shell of the supernova brightness distribution. The 
strongest knot was located northeast of the center, with 
a peak brightness temperature of ~ 3 x 10 8 K. B91 also 
showed the presence of protrusions in the brightness distri- 
bution of the supernova, at a distance of ~ 1.9 mas from the 
centre. We note that the radial distances to features A, B, 
and C in Fig. 1 are comparable to the radial distances to 
the protrusions observed by B91, and therefore one could 
be tempted to identify features A, B, and C with stationary 
clumps in the CSM of SN1986J. However, caution must be 
taken, since not only the position angles are different, but 
also our synthesized beam is significantly larger than that 
of B91. Our image also shows three spots outside the shell 
structure at contours above 3 times the rms noise off-source. 
These spots have distances and position angles with respect 
to X that are significantly different from those obtained by 
B91, and therefore do not correpond to the same physical 
features shown in B91. The peaks of these spots range from 
0.23 to 0.25 mjy/beam, less than the minimum of emission 
within the shell. Therefore, they could be mere artifacts of 
the image reconstruction procedure, implying that the pro- 
trusions seen by B91 have disappeared during the period 
ranging from September 1988 to February 1999, and point- 
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Table 1. Flux density and radial distance to X for the fea- 
tures shown in panel (b) of Fig. |l]. The position of X is a = 



02t 1 22l n 31 s 3214, S = 42° 19.'57'.'260 





Flux 


Radial distance 




mjy beam 


(mas) 


A 


1.13 


2.14 


B 


0.62 


1.10 


C 


0.77 


0.99 


D 


0.74 


1.41 



ing to a change in the density properties of the presupernova 
wind. 

Our image also shows that the radio shell morpholog 
of SN 1986J is asymmetric. The four bright features in Fig 
seem to delineate a highly distorted shell -whose minimum 
of emission is not at its centre- and are indicative of a signif- 
icant deformation of the shock front. If the point X in panel 
(b) corresponds to the centre of the explosion, we obtain the 
angular radial distances shown in Table [l] for the knots in 
the shell of the supernova. It would appear then that the su- 
pernova has expanded more along the AC axis. In addition, 
the angular radial distance XA is more than twice the dis- 
tance XC, indicative of a strong asymmetric expansion. We 
note, however, that our identification of X with the centre 
of the supernova explosion is not unique. In this respect, fu- 
ture phase-reference VLBI observations will be most useful, 
since this technique will allow a proper registration of the 
maps, thus permitting the identification of the centre of the 
explosion. 

Nevertheless, most of the results presented and dis- 
cussed below do not depend on the putative centre of 
the explosion. In particular, the brightness distribution of 
SN 1986J is clearly asymmetric, independent of the location 
of the centre of the explosion. How could such an asymmet- 
ric structure have arisen? Various models have been pro- 
posed to account for asymmetric morph ologie s in super- 
novae. Blondin, Lundqvist & Chevalier ( 1996 ) have sug- 
gested that an axisymmetric density distribution in the wind 
from a supernova progenitor leads to protrusions emerg- 
ing along the symmetry axis. These authors find that for 
a power-law supernova density profile, the flow approaches 
a self-similar state in which the protrusion length is 2-4 
times the radius of the shell, after ~ 10 yr. They cite the 
supernova remnant 41.95+575 in M82 as an example where 
such axisymmetric protrusions are likely to have emerged, 
but point out that their model is not compatible with the 
protrusions seen by B91 for SN 1986J. A case where consid- 
erable asymmetry in its over all str ucture has been observed 
is SN1987A. Gaensler et al. (1997) showed that the eastern 
and western regions of the SN1987A radio remnant were 
brighter than the northern and southern regions, and inter- 
preted this as evidence for the shocked wind being axisym- 
metric in its shape and/or density distribution. From our 
image, it follows that an axisymmetric circumstellar inter- 
action does not seem to apply for SN1986J. The protrusions 
outside the shell, on the other hand, could be formed by 
an asymmetric distributi on in the wind similar to that in- 
voked by Blondin et al. ( 1996 ), though not restricted to be 
axisymmetric. Other models that predict supernova a spher - 
ical morphologies include that of Khokhlov et al. ( 1999 ), 



which have modeled jet-induced explosions of core collapse 
supernovae. The end result is a highly aspherical supernova 
with two high- velocity jets of material moving in polar direc- 
tions, and a slower moving, highly distorted ejecta contain- 
ing most of the supernova materi al. A recent VLBI image of 
41.95+975 by McDonald et al. ( 2001 ) reveals structure re- 
sembling that of a collimat ed ou tflow, as expected from the 
model by Khokhlov et al. (1999). Although in our case one 



could identify the brightest knots A and C with a collimated 
outflow, there is also strong emission in the almost perpen- 
dicular direction BD. In addition, all four knots seem well 
confined within the supernova shell. Therefore, the existence 
of two high-velocity, well collimated jets in SN 1986J seems 
unlikely. In the case of SN1987A, a directional anisotropy 
in the supernova explosion might be indicated by the fact 
that the eastern limb of SN1987A was brighter, and was 
also increasing in brightness more rapidly, than the west- 
ern limb (Gaensler et al. 1997). Although such a directional 
anisotropy in SN1986J could explain the existence of feature 
A (which roughly corresponds to the strongest knot of emis- 
sion seen by B91), it has difficulties in explaining the change 
in position angles for the other peaks, as well as the appear- 
ance of features B and D. We therefore favour an scenario 
where the strongly asymmetric brightness distribution of the 
SN 1986J shell structure is due to the collision of the super- 
nova ejecta with a clumpy (Chugai 1993, Chugai & Belous 
1999), or filamentary wind (Weiler et al. 1990). Our image 
gives support to this model, and shows that the clumpy, or 
filamentary wind must probably be inhomogeneous to pro- 
duce a highly distorted shell. 

3.2 Deceleration of the expanding shell 

At a distance of 9.6 Mpc (Tully 1998), 1 mas corresponds 
to a linear size of ~ 1.4 x 10 17 cm ~ 0.05 pc. Based on 
8.4 GHz VLBI observations, B91 found an angular size of 
~3.7 mas (« 5.3 x 10 17 cm fa 0.17 pc) for the shell of 
SN 1986J, 5.74 yr after its explosion (assuming it took place 
on 1983.0). The corresponding mean linear velocity would 
then be ~ 14700 km s _1 for the first 5.74 yr. However, this 
velocity applies only to the protrusions found by B91, not 
to the shell. Indeed, the velocities reported in B91 were cal- 
culated for the protrusions, and assuming that these origi- 
nated in the centre at the time of the explosion. These au- 
thors also pointed out that the protrusions extended from 
the centre to twice the radius of the shell, i.e., the protru- 
sions were outside the shell. Therefore, a value of ~1.85 mas 
(fa 2.7 x 10 17 cm fa 0.09 pc) for the angular size of the 
shell of SN 1986J at epoch 1988.74 is indicated, and a mean 
linear velocity of the radio shell of ~7400kms _1 is more 
appropriate to characterize the first 5.74 yr of the expan- 
sion of SN 1986J, as has been previo usly noticed by C heva- 
lier ( |l998t ) and Houck et al. ( 199^ ). Chevalier (|l99§) also 
pointed out that such a velocity at roughly the time of the 
peak flux is evidence against the synchrotron self-absorption 
mechanism acting in SN 1986J. 

The angular size of the supernova (at the 5-rms noise 
level) along the AC direction in Figjl] is 8 ~ 4.7 ± 0.3 mas, 
equivalent to fa 6.8 x 10 17 cm fa 0.22 pc. Combining this 
angular size measurement with that obtained by Bartel et 
al. (1991) for epoch 1988.74 (9 fs 1.85 mas), we then obtain 
a mean angular expansion velocity of the shell between 29 
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September 1988 (1988.74) and 21 February 1999 (1999.14) 
of « 0.14masyr _1 , which corresponds to a linear velocity of 
~ 6300kms _1 . If we assume that SN 1986J freely expanded 
for the first 5.74 yr of its life (r oc t m ,m = 1) and then 
started to decelerate, the expansion between the two epochs 
of VLBI observations is characterized by m — 0.90 ± 0.06. 
This is a mild deceleration, and contrasts with the case of 
several other supernovae, for which a strong deceleration has 
been measured (SN1979C: Marcaide et al. 2002; SN1987A: 
Staveley-Smith et al. 1993, Gaensler et al. 1997; SN1993J: 
Marcaide et al. 1997, Bartel et al. 2000) . Thus, the angular 
size that we obtain for the shell seems consistent with the 
conclusion that the bulk of the shell structure expanded for 
the first 5.74 yr at speeds smaller than ~ 15000 km s _1 (see 
above and Chevalier 1998). 

For a standard presupernova wind velocity, v m = 10 
kms^ 1 , the linear size of SN 1986J at epoch t=16.14 yr 
implies that we are sampling the progenitor wind about 
11000 yr prior to its explosion. Since the mass loss rate of 
SN 1986J seems to have been £ 2 x 10~ 4 M Q yr -1 (Weiler, 
Panagia & Sramek 1990), the swept-up mass must have 
been M sw « 2.2Mq, and the thermal electron density ~ 
8 x 10 3 cm~ 3 (for a standard density profile of the progenitor 
wind, pes oc r~ 2 ). Since the expansion of the supernova has 
not decelerated significantly between the two VLBI obser- 
vations, it follows that the swept-up mass by the shock front 
must be much less than the mass of the ejected hydrogen- 
rich envelope, M env , as otherwise we should have observed a 
much stronger deceleration. In fact, momentum conservation 
implies that M cnv ii 12 Mq, significantly larger than M sw . 
If the hydrogen-rich mass envelope was as high as 12 Mq, 
this is a hint that the progenitor of SN 1986J was probably 
a single, massive Red Super Giant (as previously suggested 
by Weiler et al. 1990), which lost mass rapidly, but man- 
aged to keep intact most of its hydrogen-rich envelope by 
the time of explosion. This result contrasts with the cases of 
SN 1993J and SN 1979C, whose hydrogen-rich envelopes had 
masses of 0.2-0.4 M (Woosley et al. 1994, Houck & Frans- 
son 1996) and & 0.9 Mq (Marcaide et al. 2002), and whose 
progenitor stars were part of binary systems. 



3.3 Total energy and magnetic field of SN 1986J 

Since the radio emission is of synchrotron origin, we can 
estimate a minimum total energy (in magnetic fields, elec- 
trons, and heavy particles) and a minimum magnetic field 
for SN 1986J. If we assume equipartition (magnetic field 
energy is approximately equal to the total particle en- 
ergy), then the minimum total energy is (Pacholczyk 1970) 
= c w (1 + fc) 4/7 3/7 J R 9/7 L 4/7 , where L is the radio 
luminosity of the source, J? is a characteristic size, C13 is a 
slowly-dependent function of the spectral index, a, <f> is the 
fraction of the supernova's volume occupied by the mag- 
netic field and by the relativistic particles (filling factor), 
and k is the ratio of the (total) heavy particle energy to the 
electron energy. This ratio depends on the mechanism that 
generates the relativistic electrons, ranging from k ~ 1 up 
to k — ra v jm^ ~ 2000, where m p and m e are the proton 
and electron mass, respectively. 

Based on snapshot VLA observations of SN 1986J car- 
ried out on 13 June 1999 (Weiler, private communication), 



we determined an spectral index of a = —0.69 ± 0.06 
(S v oc v a ) between 1.6 and 8.5 GHz. This value is very close 
to a — —0.67 obtained by B91, and is a typical value for 
supernovae that are in the optically thin radio regime. With 
a = —0.69 and 54.9GHz = 7.2 mjy from our observations, 
we obtain a radio luminosity L„ ~ 1.38 x 10 37 (D/9.6Mpc) 
ergs -1 for the frequency range 10 7 -10 10 Hz. The value of 
the function C13 is approximately 2.39 x 10 4 (for a = —0.69, 
and v\ and V2 equal to 10 7 and 10 10 Hz, respectively). 
The filling factor, <j>, is by definition < 1. In the case of 
supernova SN 1993J, Marcaide et al. (1995a, 1995b, 1997) 
showed that the radio emitting shell of SN 1993J had a 
width, AR « 0.3 -Rout, where R ou t is the outer shell ra- 
dius, and therefore <f) « 0.66. As the characteristic size 
for SN 1986J, we will take half the largest diameter of the 
shell, 9 — 2.35 mas, which corresponds to a linear size of 
R = D ■ 9 » 3.4 x 10 17 cm. With these values, the minimum 
total energy is then 

~ 1.13 x 10 48 (1 + kf 4Q e 9l'J 5 D\T erg 

where 0o.66 = (0/0.66) 62. 35 = (0/2.35 mas) , and Dg.e — 
(D/9.6Mpc). The value of the magnetic field that yields 



E. 



B n 



(min) 



is then equal to 
lO.6(l + fc) 2/7 0, 



-2/7/J-6/7 
"0.66 (7 2.35 

.-.(tot) 



-2/7,. 
'9.6 



Since 1 < k < 2000, E^J can have values between 
~ 2 x 10 48 and ~ 9 x 10 49 erg, while the corresponding 
values of the magnetic field can lie between ~ 13 and 
~ 93 mG (for 4> = 0.66 and D = 9.6Mpc). These values 
of the magnetic field for SN 1986J are in agreement with, 
e.g., those obtained for SN 1993J at similar radii. For ex- 



ample, Perez- Torres, Alberdi, and Marcaide (2001) showed 



that B « 0.58 (r/3 x 10 16 cm)" 1 G, which results in a value 
of ~ 50 mG for a radius of 3.4 x 10 17 cm. Since it is very 
unlikely that the magnetic field energy density in the wind is 
larger than its kinetic energy density, i.e., B 2 /8ir < pv 2 J /2, 
we can then obtain an upper limit for the magnetic field in 
the circumstellar wind of SN 1986J: 



B c 



^(Mw„) 1/2 r- 1 » 0.78(M_ 4 wio) 1/2 rr 7 1 mG 



where M_ 4 = M/10" 4 M Q yr" 1 , v w = u^/lOkms" 1 , 
m = r-/10 17 cm, and we have assumed a standard wind 
density profile. For SN 1986J, M = 2 x 10~ 4 Meyr -1 and 
r = 3.4 x 10 17 cm, and we obtain B CB <, 0.32 mG. These 
arguments suggest that the magnetic field in the shell of 
SN1986J is in the range 13-93 mG, or about 40 to 300 
times the magnetic field in the circumstellar wind. Since a 
strong shock yields a fourfold increase in the particle den- 
sity, the post-shock magnetic field is also fourfold increased. 
Hence, compression alone of the circumstellar wind magnetic 
field cannot account for the large magnetic fields existing 
in SN 1986J, and other field amplification mechanisms are 
needed to be invoked, e.g., turbulent amplification (Cheva- 
lier 1982; Chevalier & Blondin 1995). The same conclusion 
was reached for the case of SN 1993J (Fransson & Bjornsson 
1998, Perez- Torres, Alberdi & Marcaide 2001), where mag- 
netic field amplification factors ~ 100 were found to be nec- 
essary. 
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4 SUMMARY 

We used 5 GHz global VLBI observations of SN 1986J taken 
on 21 February 1999 to obtain a high-resolution radio image 
of the supernova about 16 yr after its explosion (assumed to 
have occurred on 1983.0). Our main results can be summa- 
rized as follows: 

(i) The radio structure of SN 1986J is shell-like, even 
after more than 16 yr since its explosion. However, the shell 
is highly distorted, indicative of a strong deformation of the 
shock front. 

(ii) Assuming that the shell seen in our 1999.14 image 
corresponds to the shell in the 1988.74 image, the expansion 
of SN1986J has decelerated from ~7,400 km s" 1 down to 
~6,300 km s _1 . This mild deceleration can be characterized 
by a power-law (r oc t m ) with a deceleration parameter m — 
0.90 ± 0.06. Since the swept-up mass must be ~ 2.2 Mq 
(for M = 2 x 10" 4 M yr _1 and v w = 10 km s" 1 ), this 
mild deceleration seems to indicate a mass of the hydrogen- 
rich envelope ejected at explosion of <; 12 Mq. In this case, 
the supernova progenitor kept almost intact its hydrogen- 
rich envelope, in spite of its strong mass-loss wind rate, and 
is a strong hint that SN 1986 J went off in a single, massive 
star scenario. 

(iii) We show that the brightness distribution in the 
shell is strongly asymmetric, and suggest that the most 
likely scenario for its origin is the collision of the supernova 
envelope with a clumpy, or filamentary, anisotropic wind. 
However, we cannot exclude that other mechanisms, such 
as a directional anisotropy in the supernova explosion, may 
also play a role in shaping the radio emitting structure of 
SN1986J. We detect three emission spots at the 3-rms noise 
level outside the shell. Though these spots could be protru- 
sions alike those detected by Bartel et al. (1991), we suggest 
that they are mere artifacts of the image reconstruction pro- 
cedure. 

(iv) The radio flux of SN 1986J is ~ 7 mjy, which re- 
sults in a luminosity of about 1.4 x 10 37 ergs _1 . By as- 
suming equipartition between fields and particles, we esti- 
mate a minimum energy for the supernova shell in the range 
^mfn' ~ 0.2-9 x 10 49 erg, which corresponds to a minimum 
magnetic field of B m i n ps 13-93 mG. We find that the aver- 
age magnetic field in the circumstellar wind of SN 1986J is 
B cs 0.32 mG. Therefore, powerful amplification mecha- 
nisms (e.g. turbulence) should be acting in SN 1986J, since 
compression alone of the circumstellar field cannot account 
for the magnetic fields existing in the supernova shell. 

The high-resolution radio image of SN 1986J presented 
here shows that VLBI is a very powerful technique for in- 
vestigating the interaction of supernovae, because it allows 
to directly image the emission from the interaction region 
(supernova ejecta and shocked wind). Supernova SN 1986 J 
belongs to the rare group of radio supernovae that permit 
their follow-up with VLBI. Further high-resolution, sensitive 
VLBI observations of SN 1986J are necessary to provide us 
with a more detailed view of its evolution and interaction 
with the circumstellar medium. The use of phase-referencing 
will be mandatory, since this technique will not only ensure 
the detection of the supernova, but will also allow a proper 
registration of the images for all epochs. This registration 
will likely permit the unambiguous identification of the cen- 



tre of the explosion. Moreover, such observations would also 
allow to shed further light on the apparent asymmetric pre- 
supernova wind of SN1986J, obtain a better determination 
of the deceleration parameter, and definitively confirm, or 
reject, the protrusions outside the supernova shell structure 
we have marginally detected here. 



ACKNOWLEDGMENTS 

We thank the National Radio Astronomy Observatory 
(NRAO) for providing us with the archival data used in this 
paper. We thank Roger Chevalier, Nikolai Chugai, and Alan 
Roy for comments on the manuscript. We are grateful to an 
anonymous referee for many substantive and useful com- 
ments, which significantly improved our manuscript. This 
research has been supported by a Marie Curie Fellowship of 
the European Community (contract IHP-MCFI-99-1), and 
by the Spanish DGICYT grants AYA2001-2147-C02-01/02. 
KWW thanks the Office of Naval Research (ONR) for the 
6.1 funding supporting this research. The NRAO is a facility 
of the National Science Foundation operated under cooper- 
ative agreement by Associated Universities, Incorporated. 
The European VLBI Network is a joint facility of European 
and Chinese radio astronomy institutes funded by their na- 
tional research councils. 



REFERENCES 

Bartel, N., Rupen, M. P., & Shapiro, I. I. 1989, ApJ, 337, L85 
Bartel, N., et al. 1991, Nature, 350, 212 (B91) 
Bartel, N., et al. 2000, Science, 368, 610 

Beasley, A. J., & Conway, J. 1995, ASP Conference Series, Volume 
82, p. 328 Eds.: J. A. Zensus, P.J. Diamond, and P.J. Napier, 
San Francisco 

Blondin, J.M., Lundqvist, P., & Chevalier, R.A., 1996, ApJ, 474, 
257 

Chevalier, R.A., 1982, ApJ, 258, 790 

Chevalier, R.A., 1987, Nature, 329, 611 

Chevalier, R.A., 1998, ApJ, 499, 810 

Chevalier, R. A., & Blondin, J. M. 1995, ApJ, 442, 312 

Chugai, N.N., 1993, ApJ, 414, L101 

Chugai, N.N., & Belous, M.L. 1999, ARep, 43, 89 

Gaensler, B.M., et al. 1997, ApJ, 479, 845 

Houck J.C., & Fransson C, 1996, ApJ 456, 811 

Houck, J.C., et al. 1998, ApJ, 493, 431 

Khokhlov, A.M., et al. 1999, ApJ, 524, L107 

Leibundgut, B., et al. 1991, ApJ, 372, 531 

Marcaide, J.M., et al. 1995, Nature, 373, 44 

Marcaide, J.M., et al. 1995, Science, 270, 1475 

Marcaide, J.M., et al. 1997, ApJ, 270, 1475 

Marcaide, J.M., et al. 2002, A&A, 384, 408 

McDonald, A.R., et al. 2001, MNRAS, 322, 100 

Nadyozhin, D.K., 1985, Ap&SS, 112, 225 

Pacholczyk, A. G. 1970, Radio Astrophysics (San Francisco: Free- 
man) 

Perez-Torres, M.A., Alberdi, A., & Marcaide, J.M. 2001, A&A, 
374, 997 

Rupen, M.P., et al. 1987, AJ, 94, 61 
Staveley-Smith, L., et al. 1993, Nature, 366, 136 
Shepherd, M.C., Pearson, T.J., & Taylor G.B. 1995, BAAS, 26, 
987 

Tully, R.B. 1998, Nearby Galaxies Catalogue (Cambridge: Cam- 
bridge Univ. Press) 
Weiler, K. W., Panagia, N., & Sramek, R. A. 1990, ApJ, 364, 611 



© 2002 RAS, MNRAS 000, |-?? 



A distorted radio shell in the young supernova SN 1986 J 7 

Woosley S.E., et al. 1994, ApJ 429, 300 



© 2002 RAS, MNRAS 000, 



